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Q\ ■ I. INTRODUCTION 
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^ ' String theory contains moduli fields - massless fields that move on string ground state manifolds [Q . If supersym- 
metry is unbroken, these massless fields remain massless to all orders in perturbation theory, but it is assumed that all 
' ' ' moduli obtain mass through non-perturbative interactions at some high scale, with perhaps a few exceptions. Among 
CO the many moduli the dilaton is particularly interesting because it determines the string and gauge coupling. String 
' ' moduli have only nonrenormalizable couplings to light fields and their typical range of variation is the Planck scale. 
^ , Monopoles are one-dimensional topological defects that carry a magnetic charge, but their main relevant attribute 
<^ ' is that they behave as stable non-relativistic (NR) particles (see, for example, [|[||). In addition to the good old GUT 
^ monopoles, there are many stringy monopoles, dyons and other exotic creatures Since some grand symmetry is 
■ expected to brake into a lower one, monopoles and other exotics are expected to be produced via the Kibble mechanism 
' (see however [§-01). GUT type monopoles are the most dangerous objects, since they are expected to have masses of 
the order of M„ ~ 10^^ GeV, and an initial abundance of {pm/Pr)c ~ lO-^\M^/lO^^GeV){10^^GeV/T^), where 
is the monopole energy density, pr is the radiation energy density, and Tc is the temperature when monopoles were 
created 

Because of their NR nature, the energy of monopoles decreases at a slower rate than that of radiation, leading, if 
' left alone, to early monopole domination The Kibble mechanism predicts a present abundance of pm ~ lO^^pc, 
where pc is the critical energy density. But, observational limits on the presence of monopoles today, imply that the 
O fraction of monopole to the critical energy density does not exceed unity Pm/Pc < 1- This discrepancy is referred to 
O - as the monopole problem. Any other NR relics produced early enough will present the same difficulties, and since we 

' will use only NR nature of monopoles for our analysis, our results are applicable to other NR relics as well. 
^ , One class of proposed solutions argues that monopoles are produced in lower abundances or go through a phase 
' annihilation j^^-|7|. Another type of solutions is based on additional non-adiabatic expansion so that Pm/pc gets 
rS ' diluted. In order to reconcile the theoretical and observational limits, the equivalent of 27 efolds of volume expansion 
should be supplied. Inflationary models, assuming inflation does occur after monopoles were produced, can supply 
" " ■ much more than the needed 27 efolds pO|,|^,p| . 

We set out to explore the possible influence of moduli on the monopole problem. 

We show here that long periods of coherent oscillations of moduli can replace enough inflationary expansion and may 
relax the monopole problem. In general, moduli start out displaced by about a Planck distance from the global zero- 
temperature minimum of their potential, so when the universe cools down, they start to coherently roll or oscillate, 
creating particles as they do. During periods of coherent oscillations the universe is effectively matter dominated 
(MD), and the relative growth of monopole density slows down, therefore the bound on allowed initial density of the 
monopoles relaxes. If the duration of the coherent oscillation epoch is long enough, the bound is eliminated altogether. 

If moduli energy density is low and their mass lower than the Hubble expansion rate, they are essentially frozen, 
the only possible exception being the dilaton, since it develops a potential due to the existence of monopoles [pl]-p5[ . 
However, it turns out, as we will show, that the change from standard radiation dominated (RD) cosmology is small. 
The interesting situation is when moduli density is high, and deviations from standard cosmology are large. 

A similar idea has already been discusses in the context of axions, and in other work on modular cosmology, Jl6| , pTt . 

The effective equations of motion in a cosmological background with a massive dilaton included are well known. 
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2ij + 3i?2 + 1^2 _ 1^202 ^ 1 g20p ^ 

+ 3H<P + m'^(j) + ie2'^(p - 3p) = 
p + 3(H+i(/))(p + p) = 0, 



(1) 
(2) 
(3) 
(4) 



where we are using units in which nip = VIStt, where mp is the Planck mass. The conservation equation (Q) for any 
additional radiation or matter is not independent of the other equations, but we include it for completeness. Looking 
at (^, the dilaton couples to NR matter but not to radiation, since for radiation p — 3p = 0. In the case of a universe 
with radiation as the only source, (p = Const, is a solution and we retrieve the standard Friedmann-Robertson- Walker 
(FRW) cosmology. 

For completeness we also include the standard equations for any of the other moduli fields, assuming the dilaton is 
constant and fixed at the correct expectation value. 



starts to oscillate when H decreases below its mass to. 

In section II we show that for low moduli and dilaton density there are no substantial deviations from standard 
RD cosmology, in section III we analyze the case of high moduli density and in section IV we discuss our results and 
their validity. 



As already noted, the only possibly interesting field among moduli, for low density is the dilaton, because the others 
are trivially frozen. The dilaton density is given by 



We proceed to show that adding low dilaton density to a RD universe with a small amount of monopoles, does 
not substantially affect the standard evolution. We will treat each contribution to the potential separately using 
perturbation theory. 

We begin with the potential due to NR matter. For a universe with radiation and a massless dilaton, equations 




II. LOW MODULI DENSITY 



P0 = h4>^ + yi4>) 



(9) 



(0^ take the form. 



3i72 ^ 1^2 ^ lg20^^ 

4> + SH<i) = 

Pr+4:{H +l<j))pr = 0. 



(10) 
(11) 

(12) 
(13) 



The solutions of these equations are well known. 



4> — 4>o ~ Const. 

p^ = 66^200^^2 



(14) 
(15) 
(16) 



which is (up to scaling) standard RD cosmology. 
Now we will perturb the solutions. 



Pr + Sp, Sp = p,n < Pr 

H Hq + SH 
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Using (p^-p^, equations (||-||) become, to first order in the perturbations, 

QHoSH = ie^-^Vm + QHl5(j) (17) 

25H + GHqSH + 2H^S(I) = (18) 

S4, + 3Hn6^ + ie^-^Vm = 0, (19) 

We assume that all terms in the equation are of the same magnitude, and that the solutions take a power dependence 
on time, pm oc t", SH oc , S<j) oc , and we already know that Hq oc t^^. Demanding that all elements of each 
equation have the same time dependence forces a = (3 — 1 = j — 2. Since, in addition, we require Pm to describe 
monopoles in RD background, a = — |, therefore: 

The graphs in Figure 1 show that the numerical solutions of the exact equations indeed reveal such time dependence. 





FIG. 1. UP: InSH vs. Int, dashed=-l/2 sloped line. CENTER: InSp vs. Int, dashed=-3/2 sloped line. DOWN: 
In 541 vs. Int, dashed=l/2 sloped line. 
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Finally, since perturbation theory is applicable, we conclude that the presence of a massless dilaton does not 
substantially alter the evolution of monopole density in RD universe. Figure |^ shows a numerical solution for Pm/ Pr 
with and without a massless dilaton. 




FIG. 2. pm/pr vs. t: solid=with massless (f), dashed=without. 



We will now consider the regular mass term, and "shut off' the other part of the potential by putting p„i = 0. 
Equation (||) now takes the form (i) + 2iH(j) + m?(j) = 0. We are interested in the case of "slow-roll" in which the friction 



of the expansion and the potential balance each other, and acceleration is approximately zero, 
equation we obtain 



0. Solving this 



= 0(O)e^(*o-*'), K 



As long as the initial ratio of m/H is small enough, ensuring that the expansion time is shorter than the period of 
oscillations, the deviation from (/) — (I)q is small. 

Figure || shows a numerical solution for </> compared with our estimate, as well as the relative error. 
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FIG. 3. LEFT: <j) vs. t for a massive dilaton without monopoles: solid=numerics, dashed=our estimate. RIGHT: Relative 
error in our estimate. 
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The addition of mass, still keeping the low density condition, produces a small deviation from a constant dilaton, 
and therefore does not interfere with the standard evolution. 

Now we would like to consider the two potential terms together. Since the deviation of the dilaton from a constant 
in both cases is small, we have no reason to believe that the dilaton will act radically different now. The numerical 
solution shows that this is indeed a good assumption. The first graph in Figure ^ shows the relative error in estimating 
the massive low density dilaton in the presence of monopoles as a constant, the second graph shows Pmj Pr with the 
presence of a low density massive dilaton, compared with the standard evolution. 
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FIG. 4. LEFT: low density massive dilaton: in the presence of monopoles. RIGHT: ln(pm/pr) vs. Int: solid=with low 

density massive dashed= without. 

III. HIGH MODULI DENSITY 

High moduli density era is defined to commence when moduli density becomes comparable to the radiation density. 
In this era, there is no essential difference between the dilaton (assuming it is heavy enough) and any of the moduli, 
since the important element is the oscillations around the minimum of the potential rather than the coupling to 
matter. We have checked this assumption numerically in many cases. Moduli behave here as NR matter, and the 
universe is MD. 

The time dependence of the expansion of the universe with and without the moduli is different, therefore we use 
for comparison between the two not time but temperature. Using the facts that during RD H = 1.66y^— (where 
(7* is the effective number of degrees of freedom), and that in our model the evolution begins and ends with RD, we 
can use H instead of T. We are interested in how the limit on initial ratio of monopolc to radiation densities differs 
when moduli get added to the cosmic mix. 

A. A Simple Model 

The basic idea here is that since during oscillations moduli behave as NR matter, the ratio of moduli to radiation 
energy will grow roughly as the scale factor a(t), assuming that is constant or slowly varying during this period. 
The universe will very quickly reach moduli domination era in which ptot — oc a~^, and Pm/ptot oc Const. 
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FIG. 5. \n{pm/ptot) vs. \nH for various initial conditions. Solid=numerics, dashed= —1/2 sloped line. UP LEFT: 
Pm ■■ Pr ■■ P4. = I ■■ 9999 : 3. UP RIGHT: : p,. : = 1 : 9999 : 75. DOWN LEFT: p^ : p^ : p^ = 1 : 9999 : 300. 
DOWN RIGHT: p„ : p, : p^ = 1 : 9999 : 1200. 



Eventually, moduli get converted back into radiation and standard cosmology emerges, with a diluted pmj Pr- In the 
simplest model we parametrize the duration of moduli domination, and imagine the conversion into radiation as an in- 
stantaneous and completely efficient process. The graphs in figure ^ show a numerical solution of ^^{pm / Ptot) vs. In H. 
As we can see, a reasonable estimate for lalpm/ptot) is provided by a —1/2 slope for RD, and a horizontal line for 
MD. At an arbitrary point Hf, the moduli decay instantaneously into radiation and therefore once again 

In(^) ^ In(^) ~ ln(a(t)) ~ ilni ^ -^InH. (20) 
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FIG. 6. \n{p,n/ptot) vs. \nH with and without the moduli. 



The difference in the allowed initial monopole to total energy densities can be determined by the vertical distance 
between the —1/2 sloped line and the horizontal line at Hf (see Fig. A simple geometrical calculation shows that 
the relation between the duration of the moduli domination and the effective amount of volume expansion efolds it 
can replace is given by 



(Pm/ Ptot)with = S ^ {Pm/ Ptot)^i 



without ' 



(21) 



where 



a; = — In I — - 
2 \Hf 



(22) 



B. A More Realistic Model 



Now we would like to treat decay of moduli to radiation in a better way, by including a decay rate in the equations 
which will govern the duration of moduli domination period. Dimensional arguments lead to the estimate 

T = T-' ^ml/mli. (23) 

The inclusion of such a term is standard and leads to the approximate equation + iHp^ — 0, from which it is 
clearly seen that moduli really behave as NR matter. This approximate equation holds as long as we can replace (jP' 
with (</i^), which is justified when oscillations of moduli are much faster than the expansion rate. 
If we want to describe the moduli decaying into photons, or other forms of massless particles, we need to correct the 
conservation equations, moduli density decreases and the radiation density increases, 

Pr + A{H +\^)pr^Tp^, (24) 

while monopole number conservation still holds. Because the moduli's oscillations are very fast and we are using 
averaged values: 

(0) = (0) = 0. (25) 

The approximate equations are given by 
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/30 + 3Fp0 = -Tp^ 

Pr + AHpr = Tp^ 
3iJ2 ^ (i,/,2 + im2,/)2) + Ip^ « ip^ + ip^, 

where in the last equation we have assumed low monopole density. Solving for p^ we obtain 

p^ = (p^)oe-r(*-*o)(f (26) 
where to is the time when oscillations start. Solving for pr we obtain 

Pr = [{Pr)o - f rilo] (t)' + 'iTH,{^f/\ 

Therefore soon after moduli begin dominating, radiation density decreases only as R~^/'^. 

As explained, to is the time which oscillations begin, which is approximately when the moduli come into domination, 
therefore, = ^/P^i/tn'p ~ "^^/'J^p- At about trh — = T^^, most of the moduli have decayed into radiation, and 
the universe re-enters its adiabatic evolution. So between Hq = m'^/mp and Hrh — ^I'^irh = r/2, the ratio of 
monopole to radiation density decreases as R~^/^ cx H. Figure]^ shows our approximation for the behavior of Pml Pr 
with and without moduli. As before, x is the effective amount of volume efolds moduli replaces 




In H 



FIG. 7. \n{pm/pr) vs. \n{H): solid=with <j): our approximation, dashed=without <j> 

Figure ^ shows some comparisons between our approximation and an exact numerical solution for Pm/ Pr- The scale 
on the the vertical axis ln(p„i/pr) is determined by the choice of initial conditions, and H is in our standard units in 
which mr, = VTGtt. 
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C. Moduli relax the monopole bound 



Using (|28[), we want to derive some consequences regarding how moduli relax the monopole problem. 

• As was mentioned in the introduction, in order to solve the monopole problem, we need 27 efolds in volume 
expansion. If one of the moduli is to supply all 27 on its own, its mass has to be below the upper bound 

iln(|l^)>27 

< lO-^mp w 10« TeV. (29) 



• Since we do not want the monopolcs to interfere with nucleosynthesis, we should demand Trh > 10 MeV. We 
also need to keep in mind that baryogenesis has to occur later. The relationship between the moduli mass and 
the reheat temperature is: 

Hrh = 1.66Vff:^ ~ 77 = (30) 



nip 2 2mp 



which requires 
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> 10-"™p(^)i/6 « 100 TeV. (31) 

• Moduli should dilute the monopole to radiation ratio before the universe reaches monopole domination. This 
means that we need Pm/Pr at its highest point, Hi, to be less than imity 

{Pni/Pr)i = {Pml Pr)c\ jf- < 1, (32) 



where He « 10^° y^g7/T00{Tc/l0^^GeV)^GeV and {pm/Pr)c « lO-^\Mm/lO^^GeV){m^^GeV/Tc), leading to a 
lower bound on moduli mass, 

.».> 3200 (^)"'((^) GeV, (33, 

To summarize our results, here is a table with various values for moduli masses, the amount of (volume) efolds it 
replaces, its reheat temperature, and whether it acts before monopole domination, 





X 


Trh 


before md 


10"* eV 


112 


2 ■ 10"^" MeV 


no 


10" GeV 


49 


2 MeV 


yes 


10** GeV 


39 


6 • 10* MeV 


yes 


10^' GeV 


29 


2 • lO'-' MeV 


yes 


10'^ GeV 


25 


6 ■ 10"" MeV 


yes 


10'* GeV 


18 


6 ■ 10"^ MeV 


yes 



TABLE I. X values for different moduli mass. 



IV. DISCUSSION 



• The mechanism described here makes use of massive scalar fields that decay into radiation. String theory 
provides us with several candidate fields ~ moduli (including the dilaton). Instead of demanding that a single 
field provides all the 27 efolds, a possible scenario is that two or maybe more fields have different masses and 
therefore oscillate at different times, and each of them contributes a few of the needed efolds. As an example, 
consider two fields, one with mass of 10^ GeV, and the other with mass of 10^^ GeV. We will assume that for both 
fields the decay rate is given by (p3|). The 10^^ GeV mass field begins to dominate when ln{Hi) — ln(22-) ^ 11.5, 

and ends at hi{Hrh) — In(^) ^ —5.2. During this time it provides 25 efolds. Only later does the 10^ GeV 

mass field reaches domination, at ln{Hi) = —6.9 until \n{Hrh) ~ —32.9, and it provides additional 39 efolds. 
Therefore, it is possible to use use several heavier fields instead of a single lighter one. 

• The presence of moduli changes the standard adiabatic evolution of the universe, introducing a period (or 
periods) of matter domination. Looking at how many decades of temperature were indeed matter dominated 
will quantify the deviation from the standard evolution, 

H^)-H\f¥^)-hH'^), (34) 



H^)-hH^)- (35) 



and using (& 



Looking at (|2^) shows that the relation between x and the deviation is linear. We conclude that to substantially 
relax the monopole bound long periods of coherent oscillation are required. 
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We want to to compare our estimate to the numerical results. Table || shows the numerical and the estimated x 
values, and the relative error {E^ = ): 



m 


r 


^est 


Xnum 


Ex 


10"'' mp 


10"'''mp 


18.3 


9.8 


-87% 


10"** mp 


10"^*mp 


28.7 


18.6 


-54% 


10^'^mp 


10" ^"mp 


42.5 


30.3 


-40% 


IQ-^^mp 


10"'"' mp 


52.8 


39.2 


-35% 


IQ-^^mp 


10"''*mp 


63.2 


47.5 


-33% 



TABLE II. Accuracy of analytical estimates of x. 



As can be seen, our estimated values improve as the decay rate (and moduli's mass) get smaller. Also, there is 
a systematic "overshooting" (estimating too big an x). We understand this effect the following way, our calculation 
assumes that immediately as the moduli come to domination they begin to oscillate and the T term becomes effective. 
Looking at Figure ^, we see that there is a period in which moduli dominate and therefore the universe is MD 
{R oc t^^'^), yet the T term is not operative. This results in ^ cx H^^/^ instead of thus the numerical results 

give values that are lower than our estimate. 

We can try to improve our simple estimate by taking into account this rise, define Hm such that between Hi and Hm 
the universe is MD, but the moduli do not yet oscillate. Only below Hm the oscillations and the moduli decay start. 

Following [||, H„i = Such a calculation will yield a; = | = | ln(2-^^), which compared with the 

numerical values of x is consistently "undershooting" (estimating too small an x) . We have obtained an upper bound 

on X as well as a lower bound, and for the best estimated value we can take their average, Xav = ln( — ^13/12 — )■ Table 
III summarizes the relative error in the upper and lower bounds, as well as in the average, 



m 


r 




Ex (upper bound) 


Ex (lower bound) 


Ex (average) 


10"'' mp 


IQ-'^mp 


13.5 


-87% 


+11% 


-38% 


10"** mp 


10"^*mp 


21.0 


-54% 


+28% 


-13% 


10"'^mp 


IQ-^^mp 


31.0 


-40% 


+36% 


-2% 


IQ-'^mp 


10"*''mp 


38.5 


-35% 


+39% 


+2% 


IQ-^'^mp 


10"''*mp 


46.0 


-33% 


+40% 


+3% 



TABLE III. Errors in analytical estimates of x. 



To illustrate the two bounds. Figure shows the two bounds and the numerical behavior, for m = 10 ^^m. 
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FIG. 9. Approximate evolution, F 
dashed-dotted=the average. 



10 TTip solid=numerics, dotted=upper bound, dashed=lower bound, 
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